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Abstract 

The phenomenon of quantum nucleation is studied in a ferromagnet in the 
presence of a magnetic field at an arbitrary angle. We consider the mag- 
netocrystalline anisotropy with tetragonal symmetry and that with hexago- 
nal symmetry, respectively. By applying the instanton method in the spin- 
coherent-state path-integral representation, we calculate the dependence of 
the rate of quantum nucleation and the crossover temperature on the ori- 
entation and strength of the field for a thin film and for a bulk solid. Our 
results show that the rate of quantum nucleation and the crossover tempera- 
ture depend on the orientation of the external magnetic field distinctly, which 
provides a possible experimental test for quantum nucleation in nanometer- 
scale ferromagnets. 

PACS number(s): 75.45.+j, 73.40.Gk, 75.30.Gw, 75.50.Gg 







I. INTRODUCTION 



The tunneling of macroscopic object, known as Macroscopic Quantum Tunneling (MQT), 
is one of the most fascinating phenomena in condensed matter physics. In the last decade, the 
problem of quantum tunneling of magnetization in nanometer-scale magnets has attracted a 
great deal of theoretical and experimental interest.0 MQT in magnetic systems are interesting 
from a fundamental point of view as they can extend our understanding of the transition 
from quantum to classical behavior. On the other hand, these phenomena are important 
to the reliability of small magnetic units in memory devices and the designing of quantum 
computers in the future. And the measurement of magnetic MQT quantities such as the 
tunneling rates could provide independent information about microscopic parameters such 
as the magnetocrystalline anisotropies and the exchange constants. All this makes magnetic 
quantum tunneling an exciting area for theoretical research and a challenging experimental 
problem. 

The problem of quantum nucleation of a stable phase from a metastable one in ferromag- 
netic films is an interesting fundamental problem which allows direct comparison between 
theory and experiment J Consider a ferromagnetic film with its plane perpendicular to the 
easy axis determined by the magnetocrystalline anisotropy energy depending on the crystal 
symmetry. A magnetic field H is applied in a direction between perpendicular and opposite 
to the initial easy direction of the magnetization M, which favors the reversal of the magne- 
tization. The reversal occurs via the nucleation of a critical bubble, which then the nucleus 
does not collapse, but grows unrestrictedly in volume. If the temperature is sufficiently 
high, the nucleation of a bubble is a thermal overbarrier process, and the rate of thermal 
nucleation follows the Arrhenius law IV cx exp (— C//&bT), with ks being the Boltzmann 
constant and U being the height of energy barrier. In the limit of T — > 0, the nucleation 
is purely quantum-mechanical and the rate goes as Tq oc exp (—S c i/h), with S c i being the 
classical action or the WKB exponent which is independent of temperature. Because of the 
exponential dependence of the thermal rate on T, the temperature T c characterizing the 
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crossover from quantum to thermal regime can be estimated as A;bT c = HU/Sd- 

The problem of quantum nucleation was studied by Privorotskii! who estimated the 
exponent in the rate of quantum nucleation based on the dimensional analysis. Chudnovsky 
and Gunther! studied the quantum nucleation of a thin ferromagnetic film in the presence of 
an external magnetic field along the opposite direction to the easy axis at zero temperature 
by applying the instanton method in the spin-coherent-state path-integral representation. 
Ferrera and Chudnovsky extended the quantum nucleation to a finite temperature.! Kim 
studied the quantum nucleation in a thin ferromagnetic film placed in a magnetic field at 
an arbitrary angle.! 

It is noted that the previous results! - ! of quantum nucleation were obtained for ferromag- 
netic sample with the simplest form of the magnetocrystalline anisotropy energy such as the 
biaxial symmetry, and the model considered in Refs. 4 and 5 was confined to the condition 
that the magnetic field be applied along the opposite direction to the easy axis. The purpose 
of this paper is to extend the previous results of quantum nucleation in ferromagnetic sys- 
tem with simple biaxial symmetry to that of system with a more general symmetry, such as 
tetragonal and hexagonal symmetry. The generic quantum nucleation problem, however, and 
the easiest to implement in practice, is that of ferromagnets with a general structure of mag- 
netocrystalline anisotropy in a magnetic field applied at a some angle 9h to the anisotropy 
axis. This problem does not possess any symmetry and for that reason is more difficult 
mathematically. It is worth pursuing, however, because of its significance for experiments.0 
In this paper the magnetic field is applied in an arbitrary direction between perpendicular 
and opposite to the initial easy axis (z axis). Our interest in studying quantum nucleation of 
magnetic bubbles with a more general structure of magnetocrystalline anisotropy in an ar- 
bitrarily directed magnetic field is stimulated by the fact that the corresponding experiment 
would be most easy to perform and to interpret. Within the instanton approach, we present 
the numerical results for the WKB exponent in quantum nucleation of a thin ferromagnetic 
film with the magnetic field applied in a range of angles 7r/2 < 9n < 7r, where 8h is the angle 
between the initial easy axis (z axis) and the field. We also discuss the 9h dependence of the 



crossover temperature T c from purely quantum nucleation to thermally assisted processes. 
Our results show that the distinct angular dependence, together with the dependence of the 
WKB exponent on the strength of the external magnetic field, may provide an independent 
experimental test for quantum nucleation in a ferromagnetic film. Quantum nucleation (the 
description involves space-time instantons), being a field theory problem, is more difficult 
than tunneling of magnetization in single-domain particles, both at the conceptual and at 
the technical level. Therefore, this paper provides a nontrivial generalization of uniform ro- 
tation of magnetization vector (homogeneous spin tunneling) in single-domain magnetJkM^ 
to a nonuniform rotation of magnetization in bulk magnets with a more general structure 
of magnetocrystalline anisotropy in the presence of a magnetic field at an arbitrary angle. 
Compared with the tunneling in single-domain particles, a local tunneling event in a bulk 
magnet can trigger instability on a much greater scale, which leads to really macroscopic 
consequences. In experiments, it may be easier to monitor single nucleation events in a thin 
film than to detect the magnetization reversal in a nanometer-scale particle. Therefore, our 
theoretical results for a general structure of magnetocrystalline anisotropy in an arbitrarily 
directed field will be more applicable for experimental tests of quantum nucleation. Besides 
the importance from the fundamental point of view, processes of quantum nucleation and 
collapse of magnetic bubbles are potentially important for quantum limitations on the den- 
sity and long-term reliability of the data storage in magnetic memory devices and designing 
of quantum computer. 

This paper is structured in the following way. In Sec. II, we review briefly some basic 
ideas of quantum nucleation of magnetization in ferromagnets. In Sees. Ill and IV, we 
study quantum nucleation of magnetization in ferromagnets with tetragonal and hexagonal 
symmetry in an external magnetic field applied in the ZX plane with a range of angles 
7r/2 < Oh < 7r. The conclusions and discussions are presented in Sec. V. 
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II. THE PHYSICAL MODEL 



For a spin tunneling problem, the rate of magnetization reversal by quantum tunneling 
is determined by the imaginary-time transition amplitude from an initial state \i) to a final 

state \f) as 

U fi = {f\e- HT \i) = Jv{M(r,r)}exp(-S E /h), (1) 
where Se is the Euclidean action which includes the Euclidean Lagrangian density Ce as 

S E = I drd 3 rC E . (2) 



For ferromagnets at sufficiently low temperature, all the spins are locked together by the 
strong exchange interaction, and therefore only the orientation of magnetization M (r, r) 
can change but not its absolute value. For that reason the field M (r, r) is equivalent to the 
fields 9 (r, r) and (r, r), which are spherical coordinates of M. In this case the measure of 
the path integral T> {M (r, r)} in Eq. (1) is equivalent to 

JV 



J X?{#(r,r)}X?{0(r,r)} = Kmf[ (j^^j sinMMfo 



(3) 



where e = max (rfc+i — rfc) and S = M^jhq is the total spin of ferromagnet. Here 7 is the 
gyromagnetic ratio and M is the magnitude of magnetization. 

In the spin-coherent-state representation, the magnetic Lagrangian is given by 

Cjg = i Mof d^r)_\ {1 _ cog (r> r)] + E ^ 0) _ (4) 



7 \ dr 

The first term in Eq. (4) is a total imaginary-time derivative, which has no effect on the 
classical equations of motion, but it is crucial for the spin-parity effects. 0JzHl3 However, for 
the closed instanton trajectory described in this paper (as shown in the following), this time 
derivative gives a zero contribution to the path integral, and therefore can be omitted. 
The energy density in Eq. (4) is 

E(9, ( j ) ) = E a (e, ( j)) + E ex (9,(/ ) ), (5) 



where E a includes the magnetocrystalline anisotropy energy and the energy due to the 
external magnetic field, and E ex is the exchange energy 

E ex = | (d.M,) 2 = |m 2 [(W) 2 + sin 2 9 (V0) 2 ] , (6) 

where a is the exchange stiffness. .0 The magnetocrystalline anisotropy energy for tetragonal 
and hexagonal symmetry is shown in Sec. Ill and IV, respectively. In the semiclassical limit, 
the rate of quantum nucleation, with an exponential accuracy, is given by 

T cx exp [-Sg*/h] , (7) 

where Sg m is obtained along the trajectory that minimizes the Euclidean action Se- 

III. FERROMAGNETS WITH TETRAGONAL SYMMETRY 

In this section, we study the quantum nucleation of magnetization in ferromagnets with 
tetragonal symmetry in the presence of a magnetic field at arbitrary angles in the ZX plane, 
which has the following magnetocrystalline anisotropy energy 

E a (9, 0) = Kx sin 2 9 + K 2 sin 4 9 - K' 2 sin 4 9 cos (40) - M H X sin 9 cos - M H Z cos 9, (8) 

where Ki, K 2 and K' 2 are the magnetic anisotropy coefficients, and K\ > 0. In the absence 
of the magnetic field, the easy axes of this system are ±z for K\ > 0. And the field is 
applied in the ZX plane at 7r/2 < 9h < tt- Then the total energy is given by 

E [6 (r, r) , (r, r)] = K x sin 2 9 + K 2 sin 4 9 - K' 2 sin 4 9 cos (40) + -M 2 [(W) 2 + sin 2 9 (V0) 2 ] 

-M H X sin 9 cos (j)-M H z cos 9 + E , (9) 

where E Q is a constant which makes E (9, 0) zero at the initial state. By applying the 
similar method in Ref. 15, we can perform a Gaussian integration over the variable in 
the path integral and reduce the system to that with only one variable 5 (as shown in the 
following). Then it is possible to perform the rest of the calculation by using the instanton 
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method. This method simplifies the problem tremendously, compared to the problem where 
the action depended on 9 (t) and (r), though a complete mathematical equivalence to the 
initial problem is preserved. 

By introducing the dimensionless parameters as 



K 2 = K 2 /2K U K 2 = K' 2 /2K U H x = H x /H , H z = H z /H , 



(10) 



Eq. (9) can be rewritten as 

E (9, 0) = I sin 2 9 + K 2 sin 4 9 - K 2 sin 4 9 cos (40) + [(V0) 2 + sin 2 9 (V0) 2 ] 

2 AJ\ i 

—H x sin 9 cos — H z cos 9 + E , 



(11) 



where E (9, 0) = 2K\E (9, 0), and Hq = 2K\/Mq. At finite magnetic field, the plane given 
by = is the easy plane, on which E a (9, 0) reduces to 

E a (9, = 0) = i sin 2 9 + (k 2 - lf 2 ) sin 4 9-Hcos(9- 9 H ) + E . (12) 

The initial angle 9 Q is determined by [dE a (9, 0) /d6\ g _ g = 0, and the critical angle 9 C and 
the dimensionless critical field H c are determined by both [dE a {9, 0) /d9] e=e -jj =1 j = and 
[d 2 E a (9, 0) /d9 2 ] e=()c7 j =7 j c = 0, which leads to 

\- sin (20 o ) + iJsin (0 O - 9 H ) + 4 (l( 2 - 7f 2 ) sin 3 O cos0 o = 0, (13a) 
^ sin (29 c ) + # c sin (0 C - 9 H ) + 4 (Z 2 - sin 3 C cos C = 0, (13b) 
cos (20 c ) + # c cos (0 C - H ) + 4 ^2 - ^' 2 ) (3 sin 2 C cos 2 9 C - sin 4 C ) = 0. (13c) 



Assuming that 

as 

He 
sin 9 r 



K 2 -K 2 



<C 1, we obtain the critical magnetic field and the critical angle 



1 



(sin^) 2/3 + |cos^| 2/3 
1 



3/2 



4 (K 2 - K' 2 ) 

1 H — 

1 + I cot 9 H \ 2/S 



(l + |cot^| 2/3 ) 



N 1/2 



I + 1 2/3 
COt0j- 



7 H| 



(14a) 
(14b) 



In the low barrier limit, i.e., e = 1 — H / H c — > 0, by using Eqs. (13b) and (13c) we obtain 
the approximate equation for i] (= 9 C — 9 ) in the order of e 3//2 , 

"3- 



-eH c sm(9 c -9 H )+ri 2 



-H r sin 



+r]\eH c cos [9 c -9 H )-if 



- 9 H )+3{K 2 -K 2 ) sin (A9 C 
hl c cos (9 C -9 H )+A (K 2 - K'^j cos (40, 



0. (15) 



Introducing <5 = — o (|5|<^lin the small e limit), we derive the energy E (9, 0) as 



E (5, 0) = K 2 [1 - cos (40)] sin 4 (9 + 5) + H x (l- cos 0) sin (0 O + 5) 



(16) 



where E 1 (5) is a function of only 5 given by 
1- 



E 1 (S) 



-H c sin (6 C - 9 H ) + (K 2 - sin (40 c 



(<5 3 - 3c^) 



+ 
+4 



H c cos (0 C - H ) + [K 2 - K 2 ) cos (40 c ) 



(5 4 - 45 3 77 + 65V - 45 2 e) 



(x 2 - K'^j e5 2 cos (40 c ) . 



(17) 



It can be shown that in the region of n/2 < 9 H < it, < 9 C < n/2, rj and 5 are of the 
order of y/e, the second or third term in Eq. (17) is smaller than the first term in the small 
e limit. It is convenient to use dimensionless variables 



r' = e 1/4 r/r , r' = e l ' A uo Q r, 5 = 8/yfe, 



(18) 



where r 
becomes 



al\l, 



2K X 



-, and u = ^ r yK 1 /M Q . Then the Euclidean action Eq. (2) for it/ 2 < 9 H < n 



S E \5 (r', r') , (r', /)] = ^ J dr'dh' j -le 1 '* sin (0 O + y/tf) 



ds_ 

dr' 





+2K 2 sin 2 (20) sin 4 (9 + + 2H X sin 2 I ^ sin (9 + 



+ ^ 3/2 (V5) 2 + ie 1/2 sin 2 (0 O + VS) (V'0) 2 



+ 



r 3/2 K-?)}. 



(19) 



where 





cot 6h 


1/3 


1 + 


cote 


H 


2/3 



7-4 


cot 6h 


2/3 


1 + 1 


cot 0h\ 


2/3 



(20) 



In Eq. (19) we have performed the integration by part for the first term and have neglected 
the total imaginary-time derivative. In can be showed that for n/2 < Oh < 7r, only small 
values of contribute to the path integral, so that one can replace sin 2 in Eq. (19) by <fi 2 
and neglect the term including (V'0) 2 which is of the order e 2 while the other terms are of 
the order e 3//2 . Then the Gaussian integration over leads to 



J V{5(r',r')}e W (-^) , 



(21) 



where the effective action is 



S^" [S (r', r')] = hSe^rl J drW \m + \ (VS) 2 + \ - t) 



(22) 



The effect mass in Eq. (22) is found to be 



M 



\cot0 



H 



1 2/3 



2 J l+|cot6» ff 



2,/:i 



1 - e + 16K 9 + 4 



[k 2 - 



l+lcot^l 2 / 3 



+ 128^: 



2 {K 2 - <) 



2/3 



I cot e H \ 

l+|cot6» H | 2/3 



(23) 



Introducing the variables r = t 1 \J ' A/M and r = r'^/A, the effective action Eq. (22) is 
simplified as 



(24) 



For the quantum reversal of magnetization M in a small particle of volume V <C r 3 ,, M 
is uniform within the particle and 5 does not depend on the space f , Eq. (24) reduces to 



2 

S%" FM] = hSe^VMAV Jdr +\ (Vo? -t) 

The corresponding classical trajectory satisfies the equation of motion 

dH 1 m 3 2 



(25) 



(26) 



Eq. (26) has the instanton solution 



corresponding to the variation of S from 5 = at r = — oo, to 5 = V6 at t = 0, and then 
back to 5 = at r = oo. Eq. (27) agrees well with the result in Refs. 13 and 15. The 
associated classical action is found to be 

2 17/4 x 3I/4 



-hSe 5/ * 



I cot 6r l 1/6 



7 ff 



x 



l + \[K 2 -K, 



■1 \ 7-2|cot 6> H 
1+1 cat 9 H \ 



2/3 
273" 



1 - e + 16K' + A[K 2 - K' 



2 > i+|cote H | 2/3 



+ 128K 9 [K 2 -K 



|2/3 



' 2 ' i+jcot 6» H | 2/3 



(28) 



I COt 0£f 



Now we turn to the nonuniform problem. In case of a thin film of thickness h less than 
the size r /e 1//4 of the critical nucleus and its plane is perpendicular to the initial easy axis, 
we obtain the action Eq. (24) after performing the integration over the z variable, 



S e J f \5(T,T)]=hSe 3 / 4 r 2 h 




drd 2 r 



2-(I) 2+ ^ 2 +iK- j3 



(29) 



At zero temperature the classical solution of the effective action Eq. (29) has O (3) symmetry 
in two spatial plus one imaginary time dimensions. Therefore, the solution 5 is a function of 



) 1 /2 

u, where u = (p 2 + r 2 ) , and p 



1 /2 

(x 2 + y 2 ) is the normalized distance from the z axis. 



Now the effective action Eq. (29) becomes 
S e J f \5{f,r)] = A7rhSe 3 / 4 r 2 h 




j(3 

The corresponding classical trajectory satisfies the following equation of motion 



(30) 



d 2 5 2 d5 v^t 3 t2 

I 

du 2 udu 2 4 



(31) 



By applying the similar method,B! the instanton solution of Eq. (31) can be found numer- 
ically and is illustrated in Fig. 1. The maximal rotation of M is 5 max ~ 6.8499 at r = 
and p = 0. Numerical integration in Eq. (30), using this solution, gives the rate of quantum 
nucleation for a thin ferromagnetic film as 
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T Q oc exp (S E /h) 

= exp {-74.39^ 3 /V 2 /. 1 + |cOtg f / l 6 2/3 
1 Icot^r' 



3V 2 V l + |cot^| 2/3 



x 



1 - e + 16K 9 + A[K 



(k 2 - K' 2 ) 



128K 2 [K 2 - K< 



I cot e H 



2/3 



• (32) 



At high temperature, the nucleation of M is due to thermal activation, and the rate of 
nucleation follows Y T oc exp (— W min /kBT), where W min is the minimal work necessary to 
produce a nucleus capable of growing. In this case the instanton solution becomes indepen- 
dent of the imaginary-time variable f. In order to obtain W m i n , we consider the effective 
potential of the system 



U, 



eff 



d A rE 



cfr 



a 



M o 2 ((W) 2 + sin 2 £(V0) 2 )+i^,0) 



For a cylindrical bubble Eq. (33) reduces to 

1 f&X 
'o * \dp) 

From the saddle point of the functional the shape of the critical nucleus satisfies 

d 2 5 1 dS 



2 

U eff = AnK.herl dpp ± (|) + ± (Vo? - t) 



dtp 2 p dp 2 4 



(33) 



(34) 



(35) 



The solution can be found by numerical method similar to the one in Refs. 4 and 6. Fig. 2 
shows the shape of the critical bubble in thermal nucleation, and the maximal size is 3.906 
at p = 0. Using this result, the minimal work corresponding the thermal nucleation is 

W min = 41.3376fT 1 /ier 2 l . (36) 

Comparing this with Eq. (32), we obtain the approximate formula for the temperature 
characterizing the crossover from thermal to quantum nucleation as 

7-6|cot e H \ 2/3 



k B T c « 0.55 



x 



K^ 1 ' 4 |cotfl g | 1/6 
S l + |cot^| 2/3 

1 - e + 16K 2 + 4 (~K 2 



1 + - 3 (K 2 



+ | cot 9 H 



2/3 



1 



1 + | COt H 



2/3 



+ 128^2 (K 2 - K'^j - 



I cot Oh 



2/3 



1/2 



+ |cot^| 2/3 



(37) 
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To observe the quantum nucleation one needs a large crossover temperature and not too 
small a nucleation rate. Eq. (37) shows that ferromagnets with large anisotropy, i.e., small 
ration of K' 2 to K±, and small saturated magnetization are preferable for experimental study. 
In Fig. 3, we plot the Oh dependence of the crossover temperature T c for typical values of 
parameters for nanometer-scale ferromagnets: K 1 = 10 7 erg/cm 3 , K 2 = 0.1, K2—K 2 = 0.01, 
M = 500 emu/cm 3 , e = 0.01 in a wide range of angles n/2 < 9 H < n. Fig. 3 shows that 
the maximal value of T c is about 225 mK at Oh = 1.743. The maximal value of T c as well as 
Yq is expected to be observed in experiment. The similar Oh dependence of the crossover 
temperature T c was first observed in Ref. 15, while the problem considered in Ref. 15 was 
homogeneous spin tunneling in single-domain particles with uniaxial symmetry. 

IV. FERROMAGNETS WITH HEXAGONAL SYMMETRY 

In this section, we study the quantum nucleation of magnetization in nanometer-scale 
ferromagnets with hexagonal symmetry in an external magnetic field at an arbitrary angle 
in the ZX plane. Now the magnetocrystalline anisotropy energy E a (0, 0) can be written as 

E a (0, 0) = K x sin 2 + K 2 sin 4 + K 3 sin 6 - K' 3 sin 6 cos (60) 

-M H X sin cos - M H Z cos 0, (38) 

where Ki, K 2 , K 3 , and K 3 are the magnetic anisotropic coefficients. The easy axes are ±z~ 
for Ki > 0. By choosing K' 3 > 0, we take = to be the easy plane, at which the anisotropy 
energy can be written in terms of the dimensionless parameters as 

E a (0, = 0) = - sin 2 + K 2 sin 4 + (k 3 - K'^j sin 6 - H cos (0 - H ) + E , (39) 

where K 3 = K 3 /2K 1 and K 3 = K' 3 /2K 1 . 

Then the initial angle is determined by [dE a (0, 0) /d0] e = 0, and the critical 
angle C and the dimensionless critical field H c by both [dE a (0,0) /d0] e=e -jj =1 j = and 
[d 2 E a [0, 0) /d# 2 ] 0=(?c - = - = 0, which leads to 
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^ sin (20 o ) + # sin (0 O - 9 H ) + 4X 2 sin 3 O cos O + 6 (x 3 - K'^j sin 5 O cos O = 0, (40a) 

sin (29 c ) + H c sin (0 C - 6 H ) + 4X 2 sin 3 C cos 9 C + 6 3 - Z 3 ) sin 5 C cos 9 C = 0, (40b) 
cos (20 c ) + # c cos (0 C - H ) + 4Z 2 (3 sin 2 C cos 2 C - sin 4 9 C ) 

+6 (F 3 - K'^j (5 sin 4 C cos 2 9 C - sin 6 C ) = 0, (40c) 



Under the assumption that \K 2 \, 
H c and the critical angle as 



- K, 



< 1, we obtain the dimensionless critical field 



sin 9 r 



{sm9 H ) 2/3 + \cos9 H \ 2/3 
1 



-i 3/2 



1 + 



4K, 



Q[K 3 -K. 



(l + \cot9 H 



2/3^ 



1/2 



l + |cot^| 2/3 (l + lcot^l 273 )' 
I cot 9 H \ 2/3 



(41a) 



3^1 + |cot^| 2/3 
| cot 9 H \ 2/3 



+8(k 3 -K' 3 ) . 

V 7 (l + lcot^l 2 / 3 )' 

In the limit of small e = 1 — H/H c , Eq. (40a) becomes 

-eH c sin (9 C - 9 H ) + if [(3/2) # c sin (0 C - 9 H ) + 3i? 2 sin (40 c ) 
+12 3 - Z 3 ) sin 3 C cos C (5 - 8 sin 2 C ) + 77 {e# c cos (0 C - 9 H ) 
-T] 2 [(1/2) H c cos (0 C - H ) + 4K 2 cos (40 c ) 
+12 (k 3 - Z 3 ) sin 2 C (5 - 20 sin 2 9 C + 16 sin 4 C )] } = 0, 

where 77 = C — O which is small for e <C 1. By introducing a small variable 5 
(|5| <C 1 in the limit of e <S 1), the anisotropy energy becomes 



(41b) 



(42) 
0-0o 



E a (5, 0) =K 3 [1- cos (60)] sin 6 (0 O + 5) + H x (l- cos 0) sin (0 O + S) + £1 (5) , (43) 
where E\ (5) is a function of only 5 given by 



E 1 (5) 



1^ 

2 



# c sin (0 C - 9 H ) + KT 2 sin (40 c ) + a(k 3 - k'^) (5 sin 3 C cos 3 9 C - 3 sin 5 C cos C ) 
77 c cos (0 C - H ) + 7? 2 cos (40 c ) + 3 (Z 3 - K'^j sin 2 C (sin 4 C 



x (5 3 - 35 2 r]) + 



1^ 

8 



12 



-10 sin 2 C cos 2 C + 5 cos 4 C )] (<5 4 - 4<5 3 r/ + 6<SV - A5 2 e) + e5 2 [AK 2 cos (40 c ) 
+12 (Z 3 - #Q sin 2 C (sin 4 C - 10 sin 2 C cos 2 C + 5 cos 4 C ) 



(44) 



By applying the similar procedure in Sec. Ill, we obtain the transition amplitude Eqs. 
(21) and (22) by integrating out 0. For this case the effective mass is 



M = 1 + |cot0# 



2/3 



1 + -K 



^ Icot0 f - |2/3 



7 H 



X 



1 - e + 36X, 



3""^1 + Icot^l 273 
1 



K 3 -K 



I cot f -' 2/3 



7 # 



1 + Icot0# 



2/3 



+ 4K 2 1 + 120K 



(1 + 1 cot e H | 2/3 

, I cot 0^1 2/3 



G (/<, - /<Q I I -f 2 J()/v'; 



I cot H | 2/3 



i + |cot0^| 2/3 y i + |cot0 H 

-1 

1 



2/3 



l + |cot0„| 2 /V (l + |cot0 i ,| 2/3 ) : 



(45) 



and the prefactor A is 



A = 2 



I cot 



H 



,1/3 



1 + |cot0 H | 2/3 



1 2/3 



4 7- 4 lector + 2 fa _ %) 11 - 16 M fa, 



(46) 



In case of a small ferro magnet of volume V <C rjj, the result of quantum nucleation is 
Tq oc exp (—Sd/h), where the classical action for hexagonal symmetry is found to be 

017/4 x ol/4 

S cl = — HSe 5/4 \cot9 H 

5 



1/6 



2_ 7-2|cot0H| 2/3 /— — ,\ 11 - 12 loot i 
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In case of a thin film of thickness h less than the size ro/e 1 / 4 of the critical nucleus, we obtain 
the quantum nucleation as Tq oc exp (Se/K), with the classical action 
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And the crossover temperature for hexagonal symmetry is found to be 
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V. CONCLUSIONS AND DISCUSSIONS 

In summary we have investigated the quantum nucleation of magnetization in nanometer- 
scale ferromagnets in the presence of an external magnetic field at arbitrary angle. We 
consider the magnetocrystalline anisotropy with tetragonal symmetry and that with hexag- 
onal symmetry, respectively By applying the instanton method in the spin-coherent-state 
path-integral representation, we obtain the analytical formulas for quantum reversal of mag- 
netization in small magnets and the numerical formulas for quantum nucleation in thin fer- 
romagnetic film in a wide range of angles tt/2 < 6h < 7r. The temperature characterizing 
the crossover from the quantum to thermal nucleation is clearly shown for each case. Our 

14 



results show that the rate of quantum nucleation and the crossover temperature depend on 
the orientation of the external magnetic field distinctly. Therefore, both the orientation and 
the strength of the external magnetic field are the controllable parameters for the exper- 
imental test of quantum nucleation of magnetization in nanometer-scale ferromagnets. If 
the experiment is to be performed, there are three control parameters for comparison with 
theory: the angle of the external magnetic field 9h, the strength of the field in terms of e, 
and the temperature T. Our results show that ferromagnetic samples with large anisotropy 
and small saturated magnetization are suitable for experimental study the phenomenon of 
quantum nucleation. 

Recently, Wernsdorfer and co-workers have performed the switching field measurements 
on individual ferrimagnetic and insulating BaFeCoTiO nanoparticles containing about 10 5 - 
10 6 spins at very low temperatures (0.1-6K).@ They found that above 0.4K, the magneti- 
zation reversal of these particles is unambiguously described by the Neel-Brown theory of 
thermal activated rotation of the particle's moment over a well defined anisotropy energy bar- 
rier. Below 0.4K, strong deviations from this model are evidenced which are quantitatively 
in agreement with the predictions of the MQT theory without dissipationJlUll It is noted 
that the observation of quantum nucleation in ferromagnets would be extremely interesting 
as the next example, after single-domain nanoparticles, of macroscopic quantum tunnel- 
ing. The theoretical results presented here may be useful for checking the general theory 
in a wide range of systems, with more general symmetries. The experimental procedures 
on single-domain ferromagnetic nanoparticles of Barium ferrite with uniaxial symmetry^ 
may be applied to the systems with more general symmetries. Note that the inverse of the 
WKB exponent B~ Y is the magnetic viscosity S at the quantum-tunneling-dominated regime 
T < T c studied by magnetic relaxation measurements.0 Therefore, the quantum nucleation 
of magnetization should be checked at any $h by magnetic relaxation measurements. Over 
the past years a lot of experimental and theoretical works were performed on the spin tun- 
neling in molecular Mni2- and Fegli3 clusters having a collective spin state S = 10 (in 
this paper S = 10 6 ). Further experiments should focus on the level quantization of collective 
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spin states of S = 10 2 -10 4 . 

The ferromagnet is typically an insulating particle with as many as 10 3 ~ 10 6 magnetic 
moments. For the dynamical process, it is important to include the effect of the environment 
on quantum tunneling caused by phonons,00 nucleation spins,0 and Stoner excitations and 
eddy currents in metallic magnets.0 However, many studies showed that even though these 
couplings might be crucial in macroscopic quantum coherence, they are not strong enough 
to make the quantum tunneling unobservable.BllHil The theoretical calculations performed 
in this paper can be extended to the AFM bubbles, where the relevant quantity is the excess 
spin due to the small noncompensation of two sublattices. Work along this line is still in 
progress. We hope that the theoretical results presented in this paper may stimulate more 
experiments whose aim is observing quantum nucleation in nanometer-scale ferromagnets. 
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Figure Captions: 

Fig. 1 The instanton, corresponding to subbrrier bubble formation in a thin film by quan- 
tum tunneling, for r = 0, r = ±0.5, r = ±1, r = ±1.5, and r = ±2. 

Fig. 2 The shape of the critical bubble in a thermal nucleation of magnetization. 

Fig. 3 The 9 H dependence of the crossover temperature T c for rc/2 < 9 H < tt. Here, 
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Ki = 10 7 erg/cm 3 , k' 2 = 0.1, K 2 - K* 2 = 0.01, M = 500 emu/cm 3 , and e = 0.01. 
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